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Non-deterministic chaos is a form of low-dimensional dynamics which is characterized by the 
existence of a countable set of sensitive decision points (SDP's). Away from these points, the 
dynamics is well-behaved. Near these points, however, perturbations (e.g., thermal noise) may 
cause the outgoing trajectory to be chosen randomly. An example of a non-determinstic chaotic 
fT^ , system is given, and a statistical method of analyzing the resultant dynamics is developed. 
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^jQi Introduction Solutions of non-linear ordinary differential equations (ODE's) sometimes exhibit the phenomenon 

^ ' of deterministic chaos, displaying sensitive dependence on initial conditions and long-term unpredictability The 

. term "deterministic" is used because the ODE's are described by functions of the dynamical variables and time which 
are well-defined in the dynamical domain. An oft-cited examples is the kicked harmonic oscillator, which may be 

1 described by the equation 
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' where cu is the frequency of the oscillator, A is some constant, and f{x, t) is some arbitrary well-defined function. The 
solution of Eq. (|l|) is simply that of a simple harmonic oscillator (an ellipse in the phase plane) until the kick at t = nT 
occurs. At this point the solution jumps from one ellipse to another as determined by the magnitude Af{x,nT) of 
the kick at this point. For an appropriate choice of A and f{x, t), the solution x{t) of Eq. (|l|) is chaotic 
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^1^1 Let us now suppose that f{x, t) is not well-defined, but rather is a random function of time. By "random" , we mean 
■ that any point on this function has zero correlation with all other points on the function. The effect of this is to give 
the oscillator an unpredictable kick whenever t — nT. The solutions will thus be similar to those described above, 
with the important difference that the kick-induced jump from one phase plane ellipse to another is non-deterministic, 
fH and cannot be predicted. 
O ■ 

^ ' We will not examine the particular case of the randomly kicked harmonic oscillator (RKHO) here. We note, however, 
that the random function f{x,t) is expected to describable by some distribution about a mean value {f{x,t)). If this 
distribution is everywhere > 0, we are then led to the following observations: 

H ; 

5^ , 1. The dynamical trajectory after the kick at t = nT is completely uncorrelated with the trajectory before. 

2. We may make statistical predictions of what will occur after each kick based on the distribution of f{x,t). From 
this, we should be able to find some probability density that the dynamical trajectory will visit a particular 
point in phase space. 

Realization The randomly kicked harmonic oscillator, while providing a useful illustrative example, requires a 
random driving force of large amplitude, which seems rather contrived. We seek a more general paradigm where 
seemingly insignificant random fluctuations can induce the type of non-determinism seen above. To this end we 
design a "toy" dynamical system, whose orbits are all circles parameterized by their radius, r. The key feature of our 
toy system is that all of these circles share a common tangent point at the origin (see Fig. |l|). A set of ODE's whose 
solutions satisfy these criteria is easily found to be [|j 

x = y 

1 

y=7r-^^- (2) 
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We will not examine Eqs. (|^) in detail except to note that while the origin is a common point of all dynamical 
trajectories, it is not a fixed point. A particular circle of radius r (for positive x) is given by 
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[x — r) + y = r , 



(3) 



2/2 = 2xr - x^. (4) 

The equation for y may then be re-expressed as 

2xr - x'^ 1 

y = -X — r ~ x. (5) 
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In the hmit that x, y 0, we find simply that y — r. Excepting the case of the trivial circle {r ~ Q), we see that the 
origin is not a fixed point and any dynamical trajectory will reach the origin in finite time. 

Let us now add some very small random fluctuations to our system. Equations (|^) then become 
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where e <^ f and f{t) and g{t) are identically distributed random functions (we have named them differently to denote 
that they take on different functional values for the same t). Usually, if e is small enough, we ignore such fluctuations, 
as they have little effect on the dynamics. Indeed, as long as the dynamical trajectory is away from the origin, the 
average solution of Eqs. (^) will correspond to a solution of the unperturbed Eqs. (H). Near the origin, however, such 
fluctuations may have a dramatic effect, since any small change in x and y could land on a circle with any radius 
ranging from to oo. 

The dynamics of Eqs. are then expected to be similar to those of the RKHO. The solution (on average) moves 
along a circle of particular r until it reaches some neighborhood about the origin where the fluctuations become 
prevalent. The trajectory will leave this neighborhood on a circle of some different r which was randomly selected by 
the action of the fluctuations. Like the RKHO, we see that circles before and after the origin are uncorrelated, and 
that the dynamics is correspondingly unpredictable. Unlike the RKHO, the size of the fluctuations may be arbitrarily 
small, yet the resultant behavior will be equally non-deterministic (in essence, the random kick has been replaced by 
the inflnite divergence at the origin). We term this dynamical behavior non- deterministic chaos. The key feature of 
the underlying dynamical system is a point such as the origin for Eqs. (|^), which is reachable in flnite time, but is a 
common point among many dynamical trajectories. We denote such a point as a sensitive decision point (SDP). 



It should be noted here that the type of behavior described above is not unique to the particular conditions we have 
put forth. The case of noise induced instability, as described by Chen ^, also produces large-scale stochastic behavior 
in the presence of small random fluctuations. The underlying causes, however, are somewhat different. Noise induced 
instability occurs when the random perturbations drive the dynamical trajectory over the potential barrier created 
by a driving force js). Non-deterministic chaos is essentially the result of a singularity in the equations of motion (i.e., 
inflnite divergence at a point), which allows dynamical orbits comprising a flnite region of the phase plane to intersect 
at a single point. 

Statistical Properties In the case of the RKHO, the statistics of the orbit were obviously dependent on the statistics 
of the driving function. This is not so obvious in the case at hand, but we should be able to make some statistical 
predictions about the dynamics. More precisely, we would like to know the probability that a circle of given r is 
chosen when the orbit leaves some neighborhood about the origin. Let us deflne this neighborhood as a disk of radius 
(5, and note that an orbit leaving this neighborhood does so with angle 9, which we take as measured from the y-axis. 
Given that we expect the fluctuations to be isotropic, the probability density of picking a particular 9 is constant, i.e., 

p{9)d9 oc d9. (7) 

Next we note that everywhere except at the origin the Existence and Uniquenss Theorem applies to solutions of Eqs. 

, thus each circle of radius r is associated with a unique 9, and we may write 6* as a function of r. Substituting 
into p{9), we find 
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p{r)dr oc ^^dr. (8) 
or 



The probability of getting a circle between r and r + Ar is simply 



nr+Ar QQ,^ 

P{r,Ar)ozl —^dr^0{r + Ar)-9{r). (9) 



For the case at hand, we find 
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P(r, Ar) oc arccos -— — - — arccos — . (10) 

Discussion To date, there is one candidate for a physical non-deterministic chaotic system. This is a model of 
neutron star dynamics Q which exhibits a sensitive decision point and is expected to behave as described above [^,|| . 
Whether other physical examples of non-deterministic chaos exist is not clear at this time. However, it is known that 
apparently complex systems may often successfully be modeled by a low-dimensional deterministic chaotic system . 
We cautiously advance the possibility of an analogous situation for non-deterministic chaos. This is attractive, since 
the analysis and modeling of a non-deterministically chaotic signal is reduced to finding the SDP's, doing some curve 
fitting of the signal between SDP's, and performing some straightforward statistical analysis as described above. We 
are continuing to study such possibilities. 
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FIG. 1. Some examples of circular orbits of different radii, all sharing a common point at the origin. 
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